Hybrid skew scattering regime of the anomalous Hall effect in Rashba 
systems: unifying Keldysh, Boltzmann, and Kubo formalisms. 



Alexey A. Kovalev, 1 Karel Vyborny, 2 and Jairo Sinova 1 ' 2 

1 Department of Physics, Texas A&M University, College Station, TX 77843-4242, USA 
2 Institute of Physics ASCR, Cukrovarnickd 10, 162 53 Praha 6, Czech Republic 

(Dated: June 24, 2008) 

We present the analytical description of the anomalous Hall effect (AHE) in a 2DEG ferromagnet 
within the Keldysh formalism. These results unify the three linear response approaches to AHE 
and close the debate on previous discrepancies. We are able to identify a new extrinsic AHE regime 
dominated by a hybrid skew scattering mechanism. This new contribution is inversely propor- 
tional to the impurity concentration, resembling the normal skew scattering, but independent of the 
impurity-strength, resembling the side-jump mechanism. Within the Kubo formalism this regime is 
captured by higher order diagrams which, although weak, can dominate when both subbands are 
occupied; this regime can be detected by variable remote doping experiments. 
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The anomalous Hall effect (AHE) has been a subject 
of fundamental research in condensed matter physics for 
many decades. The anomalous Hall resistivity p xy de- 
scribes the transverse voltage with respect to the trans- 
port direction and depends on the spontaneous magne- 
tization M along the z direction. The origin of the 
AHE lies in the intrinsic band structure properties [H, 
and extrinsic spin-asymmetric scattering such as skew- 
scattering [3] and side-jump scattering jjjj]. 

Even though the AHE has been studied for a long time 
Q, it still remains a controversial theoretical subject due 
to the difficulty to obtain agreement between the dif- 
ferent linear response calculations within equivalent sys- 
tems 0, 0, 0, @, H, [13, U | ■ Recently, some consensus has 
been reached between the diagrammatic Kubo formalism 
[l(J and the Boltzmann approach Application of the 
Keldysh formalism to the problem is relatively new @, 0| 
and connection to the previous theories is required. Liu 
et al. |0] employ this approach but fail to reproduce the 
diagrammatic results (lTj because the gradient expansion 
of the collision integral is not taken into account [Sj. On- 
oda et al. 0] use the Keldysh technique formulated in 
a gauge invariant way; however, employment of the non- 
chiral basis representation lacks transparency and only 
allows for a numerical solution. 

In this Rapid Communication, we derive the kinetic 
equation that takes into account both the effects of the 
Berry curvature and the extrinsic effects. We solve the 
quantum Boltzmann equation analytically in the metal- 
lic (weak scattering) regime, finding the Hall current up 
to zeroth order in the impurity strength. Employing the 
chiral basis allows us to immediately identify semiclas- 
sical contributions 0] such as intrinsic, side-jump and 
skew-scattering and therefore make a systematic deriva- 
tion of the Boltzmann semiclassical approach. We also 
make a full connection to the results of the previous works 
using the Kubo formalism 0, 0, [n], 0] , hence bringing 
to an end the long standing theoretical debate within the 



weak scattering regime. In addition, we calculate the im- 
portant higher order (hybrid) skew-scattering diagrams. 
In the limit of high density and mobility, this hybrid- 
skew-scattering contribution dominates in the metallic 
regime and surprisingly has no dependence on the scat- 
tering strength but it is inversely proportional to the im- 
purity concentration. 

The method presented in the following is general, how- 
ever, in order to obtain simple analytical results that 
connect directly with other microscopic linear response 
calculations Qiji we restrict ourselves to 2D Rashba 
Hamiltonian with additional exchange field h: 



Hr = -k 2 /2m + ol-k ■ a x z — ha z + V(r) 



(1) 



where <x are Pauli matrices, tt = k — eA, A(t) = — Et de- 
scribes the external electric field and V^r) the impurities. 
Here and throughout the text we take h = c = 1. We 
employ a simplified model of impurity scattering, par- 
ticularly V(r) = VoJ2i^( r ~ r i)i wnere rj describes the 
positions of randomly distributed impurities. We also es- 
timate the spin-orbit coupling component of the disorder 
potential and show to be important only in the very high 
density regime [lil 

We start by writing the Dyson equation 
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where R, A, and K stand for retarded, advanced and 
Keldysh components of the Green's functions and self- 
energies, and the subscript labels the disorder free sys- 
tem. The symbol ® denotes a convolution (in position, 
time and spin). By considering Eg. |2l) and its conjugate, 
we arrive at the Kinetic equation [151 ]: 
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where G< /£< = (G K /f] K + G A /t A - G R /± R )/2. The 
iterative version of Eq. ([2]) corresponding to the repeated 
scattering by impurities is S = So ® [l + G ® S] , where 
Eo is the self-energy from a single scattering event in 
Keldysh space. The Keldysh component of this equation, 
gives the relation: 



l + G R (g>£ R 
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where for a single scattering event we have Eq = 0. 
Equations ([3]) and Q form a general closed set of equa- 
tions for G < and are solved in the following. 

In the presence of slowly varying perturbations, it 
is useful to perform the Wigner transformation, viz. 
the center-of-mass coordinates (X — (R,T)) and the 
Fourier transform with respect to the relative coordi- 
nates (k = (k, ui)). In such representation, the convolu- 
tion of two operators is approximated as A ® B « AB + 
| (d x Ad k B - dkAdxB^j , where we use the four vector 

notations dxdk — <9r.<9ic — dfd u and df = dr + eE<9k- 
Applying this to the Kinetic Eq. |(3]) we obtain: 
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J5± = fc 2 /2m ± A, A± = i(Gf - G A ) and G c< is G< in 
the chiral basis (r and T z are defined below; however, 
they don't affect the intrinsic current in the vanishing T 
limit). The Green's function G c< allows us to find the 
intrinsic Hall current along the x- axis: 
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where A± = ^/ (afc±) 2 + h 2 and fc 2 = 2m(ujp =p A±) de- 
scribe Fermi vectors for the lower/upper chiral bands. 

The intrinsic solution Eq. contains both the contri- 
bution at the Fermi level and from the Fermi sea, often 
referred to as a^L conductivity within the Kubo-Streda 
formalism. Our next aim is to find the contributions that 
arise due to impurity scattering at the Fermi level. We 
separate Eq. © into two parts, one is proportional to 
rip and the other is proportional to d^np. The part pro- 
portional to dujUp, i.e. the Fermi-surface, is: 



G<t A + ±<G A - G*xr 
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d k Ad x B), v = d k H , 

- np(E A E";. In deriving 

Eq. ©, one retains only the first order terms in electric 
field E and use the fact that our system is homogeneous 
and stationary (<9rG < = 0, dpG < = 0). 

To establish the connection with the several mech- 
anisms identified semiclassically when interpreting the 
AHE, we transform Eq. © into the chiral basis in which 
Hq takes the diagonal form S^HqS = lk 2 /2m — a z \, v = 
S^VqS, and: 
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where A = y/ (ak) 2 + h 2 , cos(d) — h/X and tan(0) = 
k y /k x . We first obtain the intrinsic Hall effect by disre- 
garding the collision integral in the r.h.s. of Eq. ([5j. In 
the chiral basis, only nondiagonal terms of G < give non- 
zero contributions to the intrinsic AHE and they are: 

G c <_ = ieE (iv+-d u [n F (A+ + A~)] /2+ 
(G R G R - G A G A ) v+-n F ) /2A 
-ieE {iVy+d^ [n F {A+ + A~)] /2+ 
(G R G R - G A G A ) Vy+np) /2A, 
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where A = i(G R - G A ) and f = i(E R - t A ). Note that 
dkE R ( A > = for the simple delta scatterers. We calculate 
gi?(A)< an( j Q reen ' s functions G fl( - A ' using the T-matrix 
approximation up to n^-linear terms (in the Pauli basis) 

m 
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where j = J < 

ff - S^)- 1 , with 7 = 7 r + «7 l , 0: 
calculated up to the lowest order: 
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where i/± = A± " ^+ = when w F < /i, ^ = ^~ 
when ujf > h and fco being the cutoff in the integration 
over the k vector. Note that we use the renormalizations 
u>p — > Wi? — ImT and /i — > h — lmT z in Eqs. (|9|lip which 
allows us to have purely imaginary self energies T, R< - A > . 
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while the diagonal contributions of side-jump are found 
by considering the diagram a) in Fig. Q] and by solving 
the diagonal components of the kinetic equation up to 
the second order in V . By considering the diagrams b)- 
f) in Fig. Q] and by solving the diagonal components of 
the kinetic equation up to the second order in Vb, we ob- 
tain the skew scattering contribution (the last term cor- 
responds to the diagram a) and the disorder-independent 
skew scattering [131]): 



Figure 1: Diagrams representing the averaging procedure in 
calculating E < in Eq. ifTTjl where the upper part of the plot 
corresponds to T, A and the lower part corresponds to E fl . 
The diagram a) leads to the side-jump contribution and the 
disorder- independent skew scattering [13], the diagrams b) 
and c) lead to conventional skew scattering, and the diagrams 
d), e) and f) lead to the higher order skew scattering. 
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we arrive at the the Hall conductivity |1< 



However, real parts 7 r and 7^ still appear in S < in Eq. 

CD- 

In order to find the current in Eq. ([7]) up to zeroth 
order in Vq, we transform all elements of Eqs. (|9llip 
into the chiral basis and solve the kinetic equation up to 
zeroth order in Vq. That solution is used to solve the 
diagonal components of the kinetic equation up to the 
second order in Vq. Note that the expansion of G c ^_,_ + 
starts from zero order terms in Vq (see Eq. |[6l)) while the 
expansion of G 



++/- 



starts from terms proportional to 

V ~ 2 which means that we only need to solve the non- 
diagonal components of the chiral kinetic equation up to 
zeroth order in Vb while the diagonal components of the 
kinetic equation has to be solved up to the second order. 

We find different components of G c< = G^ + G\< t + 
(7(-2) + G c s f + G c s ^ in the range (i) — h < u)p < h when 
only " — " chiral band is crossed by the Fermi level. The 
intrinsic contribution is already included in Eq. ((Sj) and 
its Fermi level part is G?< = -iEd w n F A~ («£2£<£A 

ha sin y 



4A_ 
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*r^& y ). By solving the diagonal components of the 
Kinetic equation up to zeroth order in Vq, we obtain the 
part of G c< proportional to V^ 2 : 
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where k± = (ak±) 2 + Ah 2 and cr__/ ++ = (<5o ± &z)/2- 
By solving the non-diagonal components of the Kinetic 
equation up to zeroth order in Vq , we obtain the non- 
diagonal elements of the side-jump contribution: 
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where al„ = f-(l . f „ 



and Xp = \/h 2 + 2a 2 mujp. 



In Eqs. (|13)14[) we have made a straightforward gen- 
eralization to a more general model of disorder: V(r) = 
J2i ^o^( r — r *) with Yi random and strength distributions 
satisfying <V % s = 0, ((V*) 2 ) dls = V Q 2 , {(V ^) dis = V, 3 
and ((Vo) 4 ) dis — V-t- F° r the disorder described after 
Eq. (J]), we have Vq = Vi = V% and for the white noise 
disorder we have Vi = 0. Note that this result reduces to 
the Kubo formalism result of Ref. [13] when the last term 
bracket is calculated up to zeroth order in the strength 
of the disorder. 

We repeat the same procedure in the range (ii) h < up 
when both chiral bands are partially occupied. By using 
Eq. © , within this limit we obtain that the intrinsic and 
side-jump contributions cancel each other, in agreement 
with Refs. the Fermi sea contribution vanishes 

( a xy — ) from Eq. ([6]), and the Hall conductivity is 
only non-zero for the higher order skew scattering arising 
from diagrams d)-f): 
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'xy 



>i,Y, 



TT 2 {k 2 _ 



a 2 k 2 _ k\ (fc 2 - k 2 _ )/m + 2 v /rf r +A|:(4 A_ - ki X^ 



16(a 4 + X%r/ 2 (X 2 F 



V 4 e 2 ha 2 \n\k 2 _/k 2 
"n 4 V£ Tr 2 (kl-kl) 



h 2 )/c 



(15) 

This contribution from the higher order diagrams d)-f) 
was not considered in prior calculations within the Kubo 
formalism [3, Q3 and only discussed without being cal- 
culated in Ref. |{J. We have also used the numerical pro- 
cedure of Onoda et al. 0] to verify this analytical result 
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Figure 2: (Color online.) The anomalous Hall conductivity 
ayx vs. the Fermi energy ujf and the spin-orbit coupling a 
(concentration of impurities n; = 10 11 cm -2 , exchange split- 
ting h = lOmeV and m* = 0.05m). The plot a) corresponds 
to the mobility 60m 2 /Vs or weak impurity scattering strength 
while the plot b) corresponds to the mobility 12 m 2 /Vs and 
larger impurity strength. In the inset, a\ y corresponds to the 
extrinsic skew scattering. 

that identifies this new extrinsic regime in the 2DEG with 
Rashba. Although the contribution in this regime is pro- 
portional to 1/n, it does not depend on Vq as V^" 1 , as it 
is usual for the skew scattering. The result in Eq. lfl5|) 
can be understood within the Boltzmann approach by 
writing the scattering rates in the chiral basis: 



9-JT 

,(k,k') = ^n i T*(k,k') aa .TA{k',k) 



r<5(efc 



(16) 

where f? (A) (k,k') = St(fc)(l - Vb-yW) -1 £(*;')• The 
asymmetric part with respect to k, k' in Eq. (fl6|) is 
responsible for the skew scattering and is proportional to 
UiVg A. Consequently, the Hall current should be pro- 
portional to 



skew 



where t± is the transport time for the ± chiral bands. 
The conventional skew scattering (V^ _1 order) appears 
due to the difference in the life-time for the ± chiral bands 
given by 7*. However, for the Rashba model when both 
subbands are partially occupied we have 7* = 0. In this 
limit, the asymmetry in the scattering still appears due to 
the difference in the Fermi energy renormalization for the 
± chiral bands given by 7^ and leads to a Vq independent 
contribution proportional to l/rii. 

In Fig. [21 we plot the anomalous Hall conductiv- 
ity as a function of the Fermi energy ojf and the spin- 
orbit coupling a for attractive impurities (Vq < 0). We 
take typical parameters corresponding to a high qual- 
ity 2DEG samples: the carrier concentration is in the 
range 10 11 cm -2 , the maximum spin-orbit coupling is 
5 x 10 _11 eVm and the mobilities are 12 and 60 m 2 /Vs. 
In the inset of Fig. [2] we analyze the importance of 
the extrinsic skew scattering caused by the impurity 



induced spin-orbit interaction Hso = A (<x x VV) • k 
(A = 0.052nm 2 for GaAs [ll|) that is always present in 
realistic systems. For the estimate we use the correspond- 
ing Hall conductivity [I3|, cr^ = (v-kt -v + k\). 
This conductivity becomes important for larger carrier 
concentrations and there should be a region of cross- 
over between the hybrid skew scattering and the extrin- 
sic skew scattering (some interference between the two 
effects may take place). In the limit (i) (ojf < ft), we ob- 
serve skew scattering behavior (a xy ~ l/n^Vo) when the 
inverse Born scattering amplitude r = Xjn^y^m ^> 1/ep 
(eF is the Fermi energy measured from the bottom of the 
band). For smaller r, all curves have asymptotic behavior 
reaching a sum of side-jump and intrinsic contributions 
as it can be seen from Eq. lfT4|) which represents the 
cross-over between the intrinsic-side-jump and extrinsic 
anomalous Hall effect. In the transition region to the 
limit (ii) (ojf > h), we observe a sudden drop of the 
Hall conductivity (see Fig. [2j) with a sign change. The 
hybrid skew scattering should be observable in samples 
with dopants situated closer to the 2DEG to maximize 
the impurity strength as it can be seen from the inset of 
Fig. [21 Onoda et al. analyze the region of rep ~ 1, 
finding a xy ~ a\: x scaling. This region is beyond appli- 
cability of our results which rely on the weak scattering 
limit t€f S> 1, since our approximations ignore the cor- 
rections to the conductivity ~ 1/rep and the gradient 
expansion in this regime is not fully justified. 

Summarizing, we analytically calculate the anomalous 
Hall current in a 2DEG ferromagnet with spin-orbit in- 
teraction using the Keldysh formalism. Complete agree- 
ment to the Kubo formula approach and to the Boltz- 
mann equation approach is obtained. By considering the 
higher order skew-scattering diagrams, we are able to cal- 
culate a Hall current due to a hybrid skew scattering 
mechanism which is dominant when both subbands are 
partially occupied or when the system has white noise 
disorder. This particular Hall current does not depend 
on the impurity sign and strength. 
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